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We analyze perturbative aspects of gauged matrix models, including those where clas- 
sically the gauge symmetry is partially broken. Ghost fields play a crucial role in the 
Feynman rules for these vacua. We use this formalism to elucidate the fact that non- 
perturbative aspects of Af = 1 gauge theories can be computed systematically using per- 
turbative techniques of matrix models, even if we do not possess an exact solution for the 
matrix model. As examples we show how the Seiberg-Witten solution for M = 2 gauge 
theory, the Montonen- Olive modular invariance for M = 1*, and the superpotential for 
the Leigh-Strassler deformation of M = 4 can be systematically computed in perturbation 
theory of the matrix model/gauge theory (even though in some of these cases the exact 
answer can also be obtained by summing up planar diagrams of matrix models). 
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1. Introduction 



In this paper we study perturbative aspects of matrix models as applied to non- 
perturbative dynamics of M = 1 supersymmetric gauge theories in four dimensions (ad- 
mitting a large iV description) The connection between the matrix model and the 
supersymmetric gauge theory proceeds by identifying the superpotential of the gauge the- 
ory with the potential of the matrix model. It was shown in building on previous 
work @,[5]||[7j , that the planar diagrams of the matrix model effectively compute the exact 
glueball superpotential for the associated supersymmetric gauge theory and thus yield, 
upon extremization, exact results for the gauge theory. There has been some further work 



in this direction H,@MOT2,13,14 



In some cases the planar diagrams of matrix model can be summed up exactly. This 
then gives rise to a dual geometry at the planar limit, from which one can read off non- 
trivial holomorphic information about the associated supersymmetric gauge theory. In this 
respect it is interesting to note that up to now all the cases where the supersymmetric gauge 
theory can be solved using strong/weak coupling dualities fall in the class of exactly soluble 
matrix models. In all these cases the solution takes the form of a dual geometry. However, 
in most cases (i.e. for a generic matter content and interactions) the exact solution of the 
corresponding matrix model is not available, even in the planar limit. 

But, even if the planar diagrams cannot be exactly summed, we still can resort to 
perturbative techniques of the matrix model. This yields, as noted in ||, a systematic 
instanton expansion in the gauge theory. Thus, for a large class of supersymmetric gauge 
theories for which we had no dual descriptions, we can now nevertheless compute in a 
systematic way instanton corrections to interesting holomorphic quantities. Thus, in a 
sense, we are going beyond duality, and we may hope that this will ultimately give us a 
new perspective about the meaning of duality in gauge theory and string theory. 

Perturbative techniques for matrix models are not completely trivial. This is because 
we are dealing with a gauged matrix model, and it is crucial to take this gauging into 
account properly. For vacua where the gauge symmetry is not broken, this can be easily 
taken into account by dividing by the volume of the gauge group, which simply leads to an 
overall factor. However, for vacua where the gauge group is partially broken, not only do we 
have to divide by the volume of the unbroken gauge group, we also have to deal with naive 
flat directions of the matrix fields, which are pure gauge degrees of freedom. To address 
this, we can implement the standard method of Faddeev- Popov ghosts, now applied to the 
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broken part of the gauge group. The main aim of this paper is to develop this further 
and apply it to a number of interesting examples. This will include examples where we 
know the exact solutions as well as some where we do not know how to sum up the planar 
diagrams. Since our emphasis in this paper is the applicability of perturbative techniques 
we illustrate the power of the perturbation theory, even for some of the examples where we 
do know how to sum up planar diagrams. We will consider in particular Af = 1* and Leigh- 
Strassler deformation of the Af = 4 super- Yang-Mills, as well as Af = 2 Seiberg-Witten 
geometry. 

As a byproduct of the results of this paper, which might be interesting to the matrix 
model specialists, we demonstrate how the matrix models with several eigenvalue supports 
in the large N limit can be studied by means of the planar diagram technique and estab- 
lished well-defined Feynman rules for it. (This subject is also discussed in ||15|| .) Another 
novelty which is not well explored in the matrix model literature is the possibility of filling 
not only the minima but also the maxima of the matrix potential (the "unstable" cuts), by 
virtue of the analytical continuation in the filling parameters. We demonstrate this with 
the example of the one matrix model with the cubic potential where we fill by eigenvalues 
both the minimum and the maximum. One can show that this model is equivalent to a 



particular case of the models of random paths studied in [^6[, where the solution can be 
written in terms of elliptic functions. 

The organization of this paper is as follows: In section 2 we show how gauge fixing in 
the one matrix model with the cubic potential is done, when the classical vacuum partially 
breaks the gauge symmetry. We establish the planar diagrammatic rules for this model. 
We show the importance of ghosts for matrix models in this context and relate it to the 
ghosts of the supersymmetric gauge theory. We also demonstrate that the Feynman rules 
for the multi-cut solutions have a nice geometric interpretation in terms of domain walls 
on the closed string world-sheet. In section 3 we study various examples. In appendix A 
we recall how the exact solution can be obtained in the case of the cubic superpotential 
as well as some connections with c = 1 strings on the self-dual radius. In appendix B we 
show how to setup the perturbation theory for massive vacua of Af = 1* where the rank 
of the gauge group is reduced. 

2. Gauge Fixing in Field Theory and Matrix Models 

2.1. The Problem 
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Fig. 1: The two planar two-loop diagrams, with combinatorial weight | and ^, 
that contribute to the order S 3 term in the free energy. 



To explain the setup and review the proposal of |l],[§|3[], let us start with a simple 
integral over a single M x M matrix $ 

where W(x) is a cubic polynomial with two critical points at x = a± and x = a% 

W'(x) = (x-a 1 )(x-a 2 ). (2.2) 

It was explained in |1| how to compute the genus zero free energy in this model if we put 
all the eigenvalues of the matrix $ at one critical point, say at a\. Shifting the matrix as 
$ — > ail + $ we obtain (up to a constant) 

W = ti Qa$ 2 + ^$ 3 ^ (2.3) 

with 

A = ai — <22- 

From this action we easily read off the Feynman rules: a propagator 1/A for the $ variable 
and a three-point vertex with weight 1. This gives for example the following two- loop 
contribution to the perturbative part of the genus zero free energy, with contributions ^ 
and \ horn the two planar diagrams of fig. 1 

^ = 1^ + ... (2.4) 

Here S = g s M plays the role of the 't Hooft parameter. 

According to 0, the planar limit of this matrix model can be used to obtain exact 
holomorphic quantities in the corresponding M = 1 gauge theory, which in this case is 
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simply a U(N) supersymmetric gauge theory with a single adjoint superfield and a tree- 
level superpotential trW(5>) given by (]2.2|) . For example, the effective superpotential is 
essentially given by the derivative of the J-'q(S), 

W e ff(S) = NS log(S/A 3 ) - 27rir S + N ^ y 1 (2.5) 
where the first term can be seen as coming from the contribution of the measure factor to 



the free energy JF ||17|1 . Here the variable 5" is identified with the chiral glueball field, 



S = _LtrW Q W Q . 
S2n z 

From the effective superpotential W e g one can read off non-perturbative information about 
the infra-red dynamics and vacuum structure of M = 1 theory. Thus, critical points of 
W e ff generically correspond to massive vacua in the low-energy theory. On the other hand, 
the difference AW e ff between the value of the superpotential at two different critical points 
determines the tension of the BPS domain wall separating the two vacua. 

In order to find the value of W e s at each vacuum, one should extremize it with respect 
to S and then reexpress the result in terms of the (bare) gauge coupling tq. As a result, 
one typically finds an instanton expansion, in which the n-instanton terms are fixed by 
the perturbative contributions to JF up to the n-loop order. For example, already the 
two-loop result ( [2.4] ) can be used to determine W e ff exactly up to two-instanton order. 

It is important to stress here that the rank M of the gauge group in the matrix model 
is completely unrelated to the rank A" of the gauge group in the corresponding M = 1 
theory. In order to appreciate this point, note that M enters the effective superpotential 
( P-5|) in a very complicated manner (via the S dependence), whereas the A" dependence 
is very simple (linear). In particular, the value of A" does not have to be large; the result 



(|2.5|) can be applied just as well to a U (2) gauge theory. Henceforth, we will be very careful 
to distinguish between M and N. 

Now let us proceed to a more general classical vacuum with Mi eigenvalues at a\ and 
M2 eigenvalues at 02 

*=(°i 
V a 2 

So in the matrix model we break the gauge symmetry as 

U(M) -> U{M X ) x U{M 2 ) 
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Within the string theory realization this corresponds to a background with two clusters 
of D-branes of charge M\ and M2 respectively. Taking both M\ and M2 to be large, we 
obtain a so-called two-cut solution of the matrix model. To find the perturbative expansion 
of this solution it is too naive to simply expand the matrix $ around this point. Indeed, 
if we shift 

*-(° n l M+*> (2.6) 



a 2 



and decompose the matrix $ in blocks 



*=l'J!" ! 12 | (2-7) 



$21 $22 

(where $ij corresponds to an ij string, going from the zth D-brane to the j'th D-brane) 
then the quadratic piece in the action takes the form 

I 1 

-A tr($fi + $ 2 i$i 2 - $i2$2i - $22) = g A tr ( $ ?i ~ ^22)- 

So, the kinetic terms for the "off-block diagonal" components $12 and $21 will vanish. 
This makes it problematic to keep track of the 12 and 21 degrees of freedom. 

This vanishing of the kinetic term for the off-diagonal components is not surprising 
since they are zero-modes. The original U(M) gauge symmetry still acts on the matrix 
configurations and the broken gauge transformations will transform a vacuum with two 
clusters of eigenvalues into a gauge equivalent state. More precisely, we now have a non- 
trivial vacuum manifold parametrized by the coset 

U(M)/U(M 1 ) x U(M 2 ). 

Since the action is U(M) invariant, the matrix integral will not depend on the choice of 
point on this vacuum manifold. The corresponding 2M1M2 zero-modes are exactly the 
components $12 and $21. 

The correct way to treat the semi-classical expansion, keeping track of the M\ and M2 
dependence, is by the method of Faddeev-Popov ghosts. We will see in a moment how this 
emerges both from the four-dimensional gauge theory and from the matrix model. But let 
us here remark that the role played by the ghosts is also suggested by going back to the 
topological string derivation of the matrix model as described in p]] . 

There one starts from a reduction to two dimensions of six-dimensional holomorphic 
Chern-Simons theory |I8 |. The six-dimensional open string field theory contains fields 
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of various ghost numbers that correspond geometrically to differential forms of different 
degree on the Calabi-Yau manifold. If we reduce the theory down to two dimensions, we 
find at the physical ghost level (among other fields) a gauged chiral scalar field <&(z), whose 
zero-mode is the variable <f> in the matrix integral. 

But there is also a contribution of the ghosts in this two-dimensional world-volume 
theory. One finds in particular a scalar ghost C(z) and a conjugate ghost B(z), that is a 
(1,1) form on the world- volume. Both are adjoint valued, with action 



— [d 2 z tr(BD A C + B[$,C]). 
9s J 



Since both scalars <f> and C reduce to their constant zero-modes, only the overall volume 
factor in the two-form B contributes in the path-integral. So we get an additional ghost 
contribution to the matrix integral of the form 

W ghoat = tr(B[$,C}), (2.8) 

where B, C are now anticommuting M x M matrices. Let us now explain in more detail 
the origin of this term more directly in the four-dimensional Af = 1 gauge theory and in 
the corresponding matrix model. 

2.2. Gauge Fixing in Af = 1 Supersymmetric Gauge Theory 

Consider Af = 1 gauge theory with a U(N) vector multiplet and one chiral matter 
multiplet in the adjoint representation of the gauge group. In Af = 1 superspace the field 
content of such theory is represented by a vector superfield V and an adjoint chiral scalar 
superfield <f>. Let Si nv (V, $, <f>) be the action of the superfields V and <f>, invariant under 
U (N) gauge transformations 

e v e *A e v e -*A 5 (2.9) 

where A is a chiral gauge parameter. 

Our goal will be to study (partial) gauge fixing in the functional integral 



/ 



VVV§V§ e Si=v(v,#,*) (2.10) 



by imposing a gauge fixing constraint on the adjoint scalar <f>. Implementing the standard 
Faddeev-Popov procedure, one finds: (a) that (partial) fixing of the U(N) gauge symmetry 
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leads to new anti-commuting chiral ghost superfields B and C; and (b) that the ghost action 
can be written as an F-term of the form ( |2.8| ). 

The first statement does not depend on the particular way of gauge fixing. It is 
simply related to the fact that the gauge parameter A is a chiral scalar and, therefore, the 
gauge-fixing function F = F(V, $, $) should also be a chiral superfield. namely, the gauge 
constraint should be of the form JTIJ : 

F = f, F = J (2.11) 

where / = f(x, 9) is some chiral function. As we review below, this implies that the ghost 
superfields are also chiral. 

On the other hand, the second statement above relies on the assumption that the 
gauge-fixing function F does not depend on the vector superfield V. Since, as we just 
explained, F has to be chiral we conclude that F = F(Q). In particular, a convenient choice 
of gauge is given by a linear function F(<&). Then, it follows from the gauge transformation 
of that under U(N) gauge symmetry F transforms as: 

5F = [$, A] 

Now, in order to apply the usual Faddeev-Popov method to the gauge condition ( |2.11|) , 
we introduce the functional determinant: 

A F = J VAVA S(F - f) 6(F - J) 
Inserting 1 into the path integral ( |2.10| ) in the form A^A^ 1 , we obtain 

Z = J VVV§V§ Ap 1 5(F - f) 5(F - 7) e Si ^ v '^ 

Introducing the chiral ghost fields B, C and expressing the Faddeev-Popov determinant 
A^ 1 in terms of the ghost action: 

A^ 1 = J VBVBVCVC exp tr J d 4 xd 2 9B (J^C + ^C^J 

+ tr [ d 4 xd 2 9B( 5 -^C +¥=cX 
J \SA 5A )\ 

= J VBVBVCVC exp tr J d 4 xd 2 9B{$, C] + c.c. 
VBVBVCVC e s ^ oat 



leads to the path integral 

Z = J VVV§V§VBVBVCVC e s ™+s G v+s^ OBt (2.12) 
where Sqf is the gauge-fixing action and S g h ost is given by 

^ghost = J d 4 xd 2 9 tr(fl[$, C\) +c.c. 

This is the tree- level contribution to the superpotential that we were after. Specifically, 
it shows that for a (partial) gauge fixing via imposing constraints on the adjoint chiral 
superfield $, the ghost action can indeed be written as the F-term. Moreover, the form 
of this term is exactly the same as the form of the ghost term ( |2.8|) in the matrix model 
action, which is in line with the general statement that potential in matrix model should 
be identified with classical superpotential in Af = 1 gauge theory [Q . 

2.3. Gauge Fixing in Matrix Models 

The ghost term ( |2.8|) can also be derived directly in the matrix model by gauge fixing 
the U (M) gauge symmetry that acts by conjugation on $ 

A convenient gauge choice is putting $ to diagonal form. This gives the condition 

$ y =0, % / i. 

Implementing this gauge fixing through the BRST formalism introduces exactly the above 
ghost fields; see PU[ , |2"T| ] for more discussion of ghost fields and gauge fixing in matrix 



models. 

Decomposing the ghosts also in the block form Q2.7p , we see that after the shift (|27 
the kinetic term of the ghosts is given by 

A tr(B 21 C 12 ) - A tr(B 12 C 21 ). 

So, in the case of the ghosts it is the 11 and 22 blocks that are not propagating and the 
12 and 21 block that are "physical." 

We conclude that in the reduction to the matrix integral the 11 and 22 strings represent 
physical matter fields and that the 12 and 21 strings represent ghost degrees of freedom. 



8 



This makes sense physically, since, as we already explained, in this two-cut classical vacuum 
with reduced gauge symmetry U(M\) x U(M 2 ) the matrix elements in the 11 and 22 blocks 
cannot be obtained by gauge transformations and thus they are classically not pure gauge, 
whereas the 12 and 21 blocks are pure gauge. In perturbation theory we therefore are left 
with only the ghosts in the 12 and 21 blocks. 

Before we turn to the Feynman rules that all this implies, let us point out that this 
interpretation is consistent with the multi-cut solution of the large M limit of the matrix 
integral. Here we first reduce the matrix integral to eigenvalues 

Z= fl[dX I l[(X I -Xj) 2 ew-J2 W M- ( 2 - 13 ) 

J I KJ ® s I 

In the case of a two-cut solution we can split the eigenvalues Xj in two subsets. The first 
subset of Mi eigenvalues X^ are located around the first critical point a±, the second subset 

(2) 

of M 2 eigenvalues Xj are located around the second critical point a 2 . In a semi-classical 
expansion these two critical points and the corresponding eigenvalues can be thought to 
be well-separated. We can regard the two sets {Xj } and {Xj } as eigenvalues of two 
matrices, a M\ x M\ matrix $n and a M 2 x M 2 matrix $22 with matching potentials W . 
In the saddle-point approximation after the shift ( |2.6| ) this gives the action 

w tree = tr(iA*?! + i*;o + tr(-iA$2 2 + I$y _ (214) 

From the eigenvalue representation of the matrix integral it is clear that the only way 
these matrices $n and $22 interact is through the Jacobian factor 

1, j 

(This is clearly true for arbitrary W.) This term can be exponentiated directly in the 
action (see also jyl) giving the effective action 

2 tr log ($11 ® 1 - 1 ® $22) • 

To bring out clearly the Mi and M 2 dependence, this part of the Vandermonde deter- 
minant can also be exponentiated by using the two pairs of ghosts (B 2 i, Ci 2 ) and (-B12, C 2 i). 
(We have two pairs because of the square of the Vandermonde in (|2.13| ).) In order to re- 
produce the right determinant the action of these ghosts should be 

host 

= tr(B 2 i$iiCi2 +C 21 $ 11 B 12 ) 

(2.15) 

+ tr(Si2$22C 2 i + Ci 2 $ 22 B 21 ). 
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Fig. 2: The three planar two-loop diagrams, with combinatorial weight |, 1 and 
\ respectively, that contribute to the order S1S2 term in the free energy. The grey 
propagator indicates a bosonic <&ii or $22 held; the dashed propagator indicates a 
B, C ghost of type 12 or 21. The labeling of the hole or index loop is also indicated. 

But this is exactly the action (|2.8|) restricted to the propagating fields: the 11 and 22 
blocks of $ and the 12 and 21 blocks of £?, C. 

From the two contributions to the action ( |2.14| ) and ( |2.15| ) we can read off the Feynman 
rules. We have propagators (we suppress the obvious matrix indices) 

($ii$u) = 

($22*22) = --T, 

( -B12C21 ) = ( -B21C12 ) = 

and all three-point vertices have weight 1. 

As a check of this perturbative prescription with the known properties of the two- 
cut solution we will compute in this case the two-loop contribution to the free energy 
J-"o(Si, S2). From the explicit answer to the large M solution we know this term is given 
by 0: 

i_ Q 5 3 _ bS lS 2 + bS.Si - ^ (2.16) 

The coefficients ±2/3 have already been computed. They come from the two diagrams in 
fig. 1 in which only $n and $22 (and no ghosts) propagate. 

The coefficients ±5 are given by the mixed diagrams in which also the ghosts B, C 
appear. Now there are three diagrams to consider, which are given in fig. 2. Here the 
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following factors contribute to the weight of the diagram: the symmetry factor of the 
(colored) graph, the extra minus signs of the ghost loops, the extra minus sign for the 
$22 propagator compared to the $n propagator, and the fact that their are two flavors of 
ghosts (B and C) running through each ghost loop. With these considerations taken into 
account, the three diagrams give a total combinatorial weight to the SfS^ term of 

1 • (-1) ■ 2 + 1 • (-1) • 2 + 1 • (-l)M-l) " 4 = -5. 

This indeed reproduces the second and third term in (|2.16| ). 

2.4- Relation to A2 and 0(2) models on planar graphs 

We will now argue that this two-cut model corresponds to the A2 "quiver" modelB on 



planar graphs introduced and studied in [22]. Indeed, let us consider the Feynman rules of 



the previous subsection (we choose the dimensionful parameter A = 1): if we revert at the 
same time the sign of the propagator ($22 $22) from +1 to — 1 and the sign of the weight 
of each ghost loop from +2 to —2, it is the same as to revert the sign of £2. The latter 
will lead to only positive coefficients in the formulas of the type ( |2.16|) for the expansion 
for .^(Si, S2) given in the next section. It is easy to check this statement inductively: if 
we add one ($22 $22) to any diagram (like diagrams in fig. 2) it adds up one extra loop 
weighted with the factor S2, so their sign changes are compensated. The same about a 
ghost loop: its addition leads to a new loop with the S% factor, so their sign changes are 
again compensated. 

Hence we can write down the equivalent matrix model with the potential: 



W = tr 



$2 + $ 3 + $2 + $ 3 + c t C + C tC$i + CC% 

2 3 2 3 2 



where $1 and $2 axe Mi x Mi and M2 x M2 matrices, respectively, and C = (Ci, C%) is 
a vector of two Mi x M2 rectangular complex matrix bosonic ghosts. We recognize here 
actually the A2 "quiver" matrix model with a specific matrix potential. 

In the symmetric case S = Si = —S2 this model is equivalent (only in the planar limit, 
the difference due to the uncontractible ghost loops on graphs of a nontrivial topology) to 
the 0(2) model describing the statistics of selfavoiding (ghost) loops on planar $ 3 type 
graphs, with the factor +2 for each loop (in the more general 0(n) model one has the 
weight n for each loop |T^,^3[). This model is known to describe 2D quantum gravity 
coupled to the c = 1 matter at the selfdual compactification radius. In Appendix A we 
review the full planar solution |7| of this model from the one matrix model setup. In the 
symmetric case the result is presented in terms of elliptic parametrization. 



The corresponding Coxeter diagram consists of a circle with two nodes. 
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(a) 



(b) 



Fig. 3: (a) The distribution of eigenvalues at g s = 0; (b) The dual geometry (spec- 
tral curve) at finite 't Hooft coupling. 

2.5. Multiple phases and domain walls on the world-sheet 

We would like to put the above construction into a bit more general perspective. As 
we already mentioned we are dealing with a toy model for a brane configuration where 
we have well-separated clusters of Mi, M2, . . . D-branes. In our toy matrix model we can 
see clearly how such a multi-center geometry looks like from the open and closed string 
perspectives. This might be helpful for understanding gauge/gravity dualities for these 
kind of configurations in general. 

In the matrix model at zero coupling (g s = 0) such a vacuum state is simply given by 
the distribution of the eigenvalues in groups over the critical points of W in the complex 
eigenvalue plane as sketched in fig. 3(a). The eigenvalue density is represented as a sum 
of delta-functions. 

With the use of the large N matrix model techniques we know that in the dual closed 
string picture this geometry gets modified at non-zero 't Hooft coupling . The continuous 
eigenvalue density spreads out along branch cuts in the eigenvalue plane. In this way a 
non-trivial CY geometry emerges that is essentially given by a hyperelliptic curve obtained 
as a double cover of the eigenvalue plane as sketched in fig. 3(b): 



Intuitively the following happens: if we insert a large number of eigenvalues Mj at 
the zth critical point of W this builds up a throat region in the dual geometry where the 
circumference of the neck is given by the 't Hooft coupling g s Mi. This fact that the size 



y 2 = W'(x) 2 + deformations 



(2.17) 
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Fig. 4: In a two-cut solution the 11 strings and 22 strings (here indicated in grey and 
black) will huild up world-sheet theories out of fishnet diagrams with interaction 
given by the (super)potential expanded around the relevant critical point. The 12 
and 21 strings (here indicated by dashed lines ) form self-avoiding loops, separating 
the two phases on the world-sheet. 

of the geometry is proportional to the rank of the matrix, which is a measure of the total 
number of degrees of freedom, should be thought of as a version of the Bekenstein-Hawking 
geometric entropy, and it would be interesting to develop this interpretation further. 

We have seen that in the open string picture the character of the ij strings, stretching 
from the zth to the jth D-brane, is very different depending on whether j = i or j ^ i. The 
diagonal ii strings have interactions among themselves that are given by the expansion of 
the superpotential W around the ith critical point and can therefore be of arbitrary order. 
These interactions build up the fishnet double-line Feynman diagrams that in the large N 
limit will describe the closed string world-sheet propagating in the local geometry around 
the zth D-brane, just as in the case of a single center geometry. 

The interactions of the off-diagonal ij strings with j '< ^ i do not depend on the potential 
W(<fr). They are given entirely by the cubic interaction ( |2.15| ) that is dictated by gauge 
invariance. Note that the action is quadratic in these % ^ j strings — ghost number is 
conserved — and therefore the ghost loops will form well-defined demarcation lines on 
the closed string world-sheet separating the "phase" where the string is propagating in 
the background of the ith D-brane from the phase where the string propagates in the 
background of the jth D-brane, as sketched in fig. 4. Because the absence of interactions 
among the ij strings these loops are self-avoiding. 
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In this way we observe that the multi-cut solutions of the matrix model translated 
into a closed string picture naturally describe a system of dynamical domain walls on the 
world-sheet. These domain walls connect different conformal field theories as was analyzed 
in ||24j|. In the open string channel the domain wall corresponds to an ij string stretching 



from one throat to another. This picture of different phases of the world-sheet of a single 
closed topological string is a further application of the ideas in [IT7 . 



3. Examples 

In this section we illustrate how matrix perturbation theory can be used to obtain 
non-perturbative instanton effects in various supersymmetric gauge theories. We start 
with some familiar examples, which include Af = 2 Seiberg-Witten theory and M = 1* 
theory, where the exact answer is known to all orders. Despite the existence of the exact 
solution in these models, we will not need it here. Instead, our goal is to reproduce it by 
computing simple planar diagrams in the corresponding matrix model. 

Of course, the real power of the perturbative technique is in those models where exact 
solution is not available. It is easy to come up with simple examples of such models. A 
particular example that we discuss in this section is a massive deformation of the Leigh- 
Strassler theory, which in turn is an (exactly marginal) deformation of the M = 4 super- 
Yang-Mills |[25|| . The case that we consider corresponds to a simple 3-matrix model with 
cubic interactions, solution to which is not known even in the planar limit. Nevertheless, 
one can systematically obtain instanton corrections to the effective superpotential from 
matrix perturbation theory. Similar perturbative analysis can be applied essentially to 
any M = 1 theory that admits a large N limit. 

3.1. Seiberg-Witten Solution from Multi-Cut Matrix Models 

The fact that one can obtain the Seiberg-Witten solution from a perturbative analysis 
of the gauge theory, which in turn gets reduced to planar computations of a matrix model 
has already been noted in [R] as an interpretation of the string inspired derivation of 



Seiberg-Witten geometry in [26]. Our aim in this section is to show that even if the 
exact solution of matrix model were not available we could have nevertheless obtained a 
systematic instanton expansion for quantities of interest. So in this section we are tying 
one hand behind our back. 
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The basic idea of |26[] is to deform M = 2 theory to M = 1 by a polynomial tree-level 
superpotential W(x), which freezes the eigenvalues of the adjoint field $ to a particular 
point on the Coulomb branch. For example, in the case of U(2) gauge theory one deforms 
by a cubic superpotential of the form ( |2.2|) : 

W/reeO) = e(x - a)(x + a). 

Here we explicitly introduced the deformation parameter e, such that e = corresponds 
to the undeformed M = 2 theory. Choosing the configuration where one eigenvalue of $ 
is at +a and the other is at —a determines a point on the Coulomb branch of the original 
Af = 2 theory, and breaks the gauge group to an abelian subgroup, 

U{2) -> U(l) x U(l). 

This leads us precisely to the situation discussed in the previous section, where we 
studied vacua of M = 1 field theories with (partial) gauge symmetry breaking. Therefore, 
one should be able to compute all holomorphic quantities from the genus zero free energy 
FoiSi, S 2 ) of the corresponding two-cut matrix model. Evaluating the two-loop Feynman 
diagrams in the previous section we found the leading perturbative behaviour of the genus 
zero free energy in the two-cut matrix model with a cubic interaction: 

JT rt (Si, S 2 ) = ± Qs? - 5SfS 2 + 5S 1 Sl - ~Sf) + . . . 

One can go further and systematically compute higher-order corrections. In this way one 
finds a series expansion 

MSu S 2 ) = -\ J2 Sf log (^j + (S 1 + S 2 f log (£j + 

+ ^ (f s? - s fsfs 2 + ^p^f - ^p^f + - 5 pi) + ■ • ■ 

Here the first term receives a contribution from the measure of the unbroken gauge group 
U(Mi) x U(M 2 ) [pTTH , where each factor gives a standard term Sf log Si that reproduces 
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the Veneziano-Yankielowicz superpotential. The one-loop diagrams for $ and the ghosts 
B, C account for the A dependence of the first two terms in ( |3.1| ) 

]-S 2 -2S 1 S 2 + ~Si) log A. 

Finally, the A dependence reflects the ambiguity in the cut-off of the full U (Mi +M2) gauge 
group and should therefore multiply (Si + S2) 2 ■ The higher order perturbative terms have 
the combinatorial meaning we explained in the previous section. For example, the terms 
that involve only S\ or S2 enumerate planar cubic diagrams and were computed in |27[ . 



Note that the function J-o(Si, Sz) is symmetric in Si and — S 2 - This reflects the 
symmetry of the potential: we can exchange the stable and unstable critical points if we 
change the overall sign of the potential by g s — > —g s . Since Si = g s Mi this gives Si <-> — S 2 - 
From the combinatorial point of view this was explained in section 2.4 in terms with the 
connection to the 0(2) model on a random surface — it is an obvious property of the 
Feynman rules. 

We should now extremize the effective glueball superpotential 

W eS (S) = { Ni ^S^ " 2mT ° S ) ( 3 - 2 ) 

i 1 

In the present case we have N\ = N 2 = 1 and we will also set to zero the bare coupling tq. 

The physical quantity to compute in this model is the matrix of the U(l) x U(l) 
couplings in the effective low-energy theory. It is given by the second derivatives of matrix 
model free energy 

_ d 2 f (S) 

TlJ ~ ~ds~ds~- (3 - 3) 

Note that by a scaling argument the matrix does not depend on the deformation 
parameter e and therefore it should reproduce the coupling constant of the Af = 2 Seiberg- 
Witten theory at the relevant point of the Coulomb branch. Minimizing the effective 
superpotential ( |3.2|) , that in this case simplifies to 



gives the condition 
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So we see that at the extremum r,j takes the form 



Tn T12 \ _ f 1 -1 
T21 T 2 2 / V — 1 1 



(3.4) 



where r is the effective gauge coupling for the 'off-diagonal' U(l) C SU(2) C f/(2). Note 
that we automatically managed to get rid of the diagonal U(l) factor by setting the bare 
coupling constant to zero in eq. (|3.2|). 

The extremization of W e s (S) we can do using the perturbative expansion of To ( |3.f|) . 
However, before we do this, let us recall that, in terms of the exact solutions, this extrem- 
ization has a clear geometric interpretation 0,f26[. The free energy To can be described 



in terms of the dual geometry (|2.f 7|) that in this case of a cubic superpotential takes the 
form of a genus one curve 

y 2 = (x 2 - a 2 ) 2 + b x x + bo. (3.5) 

Here the coefficients bi, bo are determined by the 't Hooft couplings Si, 62- In particular 
one has the simple relation b\ = — 4(Si + S2). Minimizing W e ff (S) with respect to Si and 
£2 gives the condition 

St = -S 2 . (3.6) 



Therefore the algebraic curve (|3.5| ) reduces to nothing but the Seiberg-Witten curve for 
SU{2) theory [f2§: 

y 2 = (x 2 - u) 2 + A 4 , 
where one has to make the identification of parameters (with A = 2a) 

u= 1 -(tv$ 2 ) = 1 -A 2 . (3.7) 

So at the extremum the free energy To can be thought of as a function of only one variable 
S = Si = —S2 that is determined by the parameter A (or u) of the SW curve. 

Both in the matrix model and in the SW solution the exact expression for the cou- 
pling r of the off-diagonal U(l) that appears in (|3.4| ) follows directly from this geometric 



interpretation as the modulus of an elliptic curve. Given the parametrization of this curve, 
we can expand r in terms of the variable u or A and obtain the exact result 



A\ /A\ 4 /A\ 8 62048 /A x 12 



r( W ) = 21og(-j+20^-j +538(^-) + ( - ) +... (3.S) 
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We can now reconstruct this exact solution in perturbation theory by simply evaluat- 
ing the second derivative of JF at the critical point up to a fixed number of loops. Given 



the perturbative expansion ( |3.1| ) of JF in terms of a loop expansion of planar diagrams, 



we should first compute S = S± = —S2 at the extremum. This gives a series of the form 
S /A\ 4 /A\ 8 . /A\ 12 ... /A Nlf> 



A3 V a; +6 {a) +uo {a) +4620 UJ (a9) 

Plugging this into 9 2 JF / dS 2 gives us a systematic approximation of the effective coupling 
r. It is instructive to see how the instanton expansion of r computed from the n-loop free 
energy of the matrix model for various n gives a sequence of series expansions gradually 
converging to the exact result: 



A 

A, 

A\ _ /A\ 4 „ _/A\ 8 / A N 



Tl-loop = 2 log 



r 2 -ioo P = 21og(-j+20^-j +120^j + 1080 ^ - J +... 

A\ M\ 4 /A\ 8 /A xl2 



T 3 -ioop = 21og(-j+20^-j + 538 ^-j +7816(^-j +... 

A\ „„M\ 4 __^A\ 8 62048 /A x 12 



(3.10) 



r 4 - lo op = 21og(-l+20l-l +538 ( - J + — ( - ) !. 



A\ /A\ 4 /A\ 8 62048 /A x 12 



W t = 21og^-j + 20^-j + 538 ^-j +—{£) + 

As an aside we point out that the condition Si = —S2, that naturally emerges from 
minimizing the effective superpotential, means that from the point of view of the matrix 
model we are dealing with a symmetric filling of the two cuts. The exact solution to this 
model has interesting properties and is further analyzed in Appendix A. In particular there 
it is discussed that this model, as well as its generalisation with asymmetric filling of the 
two cuts, has a non-trivial scaling limit in the universality class of the c = 1 string. 

Remembering the relation Si = g s Mi, we see that because of the minus sign in (|3.6Q 
in the symmetric filling the number of eigenvalues in the unstable cut (the maximum of 
the potential) is negative. This is clearly an unphysical solution and should be interpreted 
as obtained by analytic continuation. In fact, if we put a positive number of eigenvalues 
at an unstable critical point the eigenvalue cut will not lie on the real axis but the cut will 



18 



rotate itself along the imaginary axis. (This can be seen by simply analytically continuing 
$ —>■ i$ in the Gaussian approximation.) Instead of working with negative numbers it is 
perhaps better to think of this solutions in terms of "eigenvalue holes" obtained by filling 
the Dyson sea almost to the top of the potential. 

Finally, let us point out that using matrix model results we could also obtain other 
holomorphic quantities, such as the SW periods a and a£>. For example, for the expectation 
values (tr$ fc ) one finds a nice expression: 



written in terms of the 1-form h = W"(x)dx/y, which can be interpreted as the smeared 
density of the eigenvalues of the adjoint field $ In particular, the case k = 1 gives 

rise to the SW period a. 

3.2. Af=l* Theory 

As another illustration of the perturbative technique in the matrix model applied to 
non-perturbative gauge theory, we consider a massive deformation of M = 4 gauge theory, 
the so-called Af = 1* theory. In Af = 1 superspace this theory is described by U(N) gauge 
theory with three adjoint chiral superfields and a tree-level superpotential: 



For simplicity, we also assume that the eigenvalues of all the Higgs fields are in the same 
classical vacuum (perturbation theory around other vacua is discussed in Appendix B). 
Computing planar Feynman diagrams up to 3 loops in the corresponding matrix model 
we will be able to reproduce the leading terms in the (exact) effective superpotential of 
Af = 1* theory. 

From the topology of planar Feynman diagrams in this matrix model it is easy to see 
that the free energy, !Fo(S), has the following structure 





(3.11) 




(3.12) 



k 



Following the notations of 0, henceforth we set g = 1. 
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Fig. 5: Two-loop contribution to the free energy in the 3-matrix model corre- 
sponding to the M =1* theory. The numbers next to propagators label the choice 
of one of the three matrix fields. 



Given the matrix model free energy !Fq(S), one can compute the effective superpoten- 
tial W e ff(S) using the relation (|2.5|). Furthermore, integrating out the field S in W e s(S) 
gives the effective superpotential as a function of the coupling constants. For the Af = 1* 
theory the answer can be computed explicitly H by using the matrix model techniques 



developed in [p9[. Specifically, one obtains 



W eS = --^ r E 2 (r), (3.13) 



where r = tq/N and E2(r) is the Eisenstein series. This agrees with the analysis of 
based on field theory dualities. Up to an additive constant, we can write the effective 
superpotential ( p.l3| ) as a power series in the variable q = exp(27rzr): 



= Nm s (q + 3q 2 + 4q s + 7q 4 + 6q 5 + ...), (3.14) 



Our goal is to reproduce this result by the perturbative technique in the corresponding 
3-matrix model 

r 3 

/ d$ exp-tr($ 1 [$ 2 ,$ 3 ] + y^^ 2 ). (3.15) 

J i=i 

Namely, computing the planar Feynman diagrams up to three loops we shall find numerical 
coefficients in the perturbative series (|3.12|) and, in particular, to check the first few 
coefficients in eq. (|3.14|) . 

The two-loop contribution to JF comes from the Feynman diagrams of the type shown 
in fig. 5. It is one of the diagrams that appears in a simple 1-matrix model with cubic 
potential, see fig. 1. The second type of 2-loop diagrams in fig. 1 does not appear here due 
to the index structure of the cubic interaction. Thus, we obtain the two-loop coefficient 
c 2 = -1. 
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2 



Fig. 6: Two types of 3-loop diagrams that contribute to J-q with one of the possible 
labeling of the propagators. 



At the next order, i.e. at three loops, there are two types of diagrams which are 
presented in fig. 6. Taking into account also the index structure of the diagrams one finds 
many different terms. Adding all of them together gives cs = 7/2. 

Summarising, up to three loops the perturbative expansion of J-'o(S) has the form: 

^o{3) = ~ + ~ + ... (3.16) 
ui 2 m° 

Note, that the expansion we find indeed has the general structure expected in ( |3.12| ). The 
relative minus signs in this expansion are due to the interaction vertices with both positive 
and negative weight arising from the commutator in ( |3.15|) . 

Substituting ( |3.16|) into fl2.5|) we obtain the leading behavior of the effective superpo- 
tential 

W eS = NSlog(S/m 3 ) - 2mr S - 3N— + UN— + ... (3.17) 

Now, integrating out the gluino field S we obtain the final expression for the effective 
superpotential 

W eS = -Nm 3 q - 3Nm 3 q 2 - 4Nm 3 q 3 + ... (3.18) 

The leading coefficients in this expression agree with the first coefficients in the expansion 
of the exact answer ( |3.14| ) , written in terms of the Eisenstein series Ei (r ) . 

Since we can do this calculation order by order, and since n-loop diagrams give rise to 
n-instanton terms in W e g, it is instructive to look at the higher order terms and to see how 
the result depends on n. For example, if we keep only the leading term in the perturbative 



series Fq, the superpotential (|3.17|) looks like: 

Weff = ^log(^)-27rzr ^-3^ 2 (3.19) 
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This leads to the effective superpotential 

Wes = Nm 3 (q + 3q 2 + 18q 3 + UAq 4 + 1350q 5 + . . .) (3.20) 

where we retained the terms of higher order in g, most of which can not be trusted in this 
approximation. 

If we compute perturbative free energy Tq to three loops, as we did above, we obtain 
the effective superpotential ( |3.18] ), where one can trust three leading terms. Moreover, the 
values of the higher order terms in ( |3.18|) are slightly "improved" compared to (|3.20|) . One 



(3.21) 



can continue and do a similar calculation up to four loops and so on. As a result, one finds 
a sequence of instanton expansions which gradually approach the exact answer (|3.14| ): 

Wi_i oop = Nm 3 q 

W 2 _ loop = Nm 3 ( q + 3q 2 + 18q 3 + UAq 4 + 1350^ + ^^q 6 + + 
V 5 5 

WVioop = Nm 3 (q + 3q 2 + Aq 3 - lQ8q 4 - 1548? 5 - + ™™ q 7 + . 

V 5 5 

W 4 _ loop = Nm 3 ( q + 3q 2 + Aq 3 + 7q 4 + 1212^ + 1^1^ + 874744^ + 

at 3/ o 2 , 3 » 4 n 5 72516 6 1657856 7 \ 
W 5 _i oop = Nm 3 ( q + 3q 2 + Aq 3 + 7q 4 + 6g 5 —q 6 q 7 + . . . J 

WV loop = Nm 3 ( q + 3q 2 + Aq 3 + 7g 4 + 6g 5 + 12g 6 + 190976 ? 7 + . . . ) 

^exact = Nm 3 ( q + 3q 2 + Aq 3 + 7q 4 + 6q 5 + 12g 6 + 8q 7 + . . . ) 

Here, the underlined terms represent the exact terms in the instanton expansion whose 
coefficients "stabilize" beyond a certain order. It is curious to note, that although all the 
numerical coefficients in the exact superpotential W exac t are integer numbers, it is not the 
case for the result obtained from a finite number of loops in matrix perturbation theory. 
Moreover, the n-loop approximation to W e xact is n °t a modular form, and one can see from 
the examples listed above that in the truncation to n loops the mistake in the (n + l)th 
coefficient is quite large. This emphasizes the fact that the Montonen-Olive duality is not 
put in by hand in this formalism, but rather is derived. In this sense, we are going one 
step beyond duality. 

Note that we can express the S'-duality of the N = 1* theory succinctly as the state- 
ment that the effective glueball superpotential W e ff(S) is given by a Legendre transform 
of a modular form, in this given by £ ! 2( r ) (with r = tq/N). 
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Fig. 7: The Feynman rules in the perturbative 3-matrix model corresponding to 
the massive deformation of the Leigh-Strassler model. 



3. 3. Massive Deformation of the Leigh-Strassler Model 

So far we considered only examples for which exact solution was already known. This 
was helpful for establishing some confidence in the perturbative technique since it did not 
rely on the existence of the exact results, which we used only to verify the perturbative 
answer. As we explained in the introduction, in most of the models we don't have this 
luxury and, therefore, perturbative analysis remains as the only tool for obtaining non- 
perturbative results, such as instanton expansion of the effective superpotential. Here, we 
consider one such model. 



Specifically, we consider a Leigh-Strassler deformation J25| of the model discussed in 
the previous subsection: 

i=l i=l 

The corresponding 3-matrix model with action given by Wtree(^i) can be solved in the 
large M limit if either g = or h = 0, but the exact solution is not known when both 
deformation parameters, g and h, are non-zero. On the other hand, perturbation theory 
is very simple, with the Feynman rules summarized in fig. 7. 
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At the 2-loop order, we find the following expression for the genus zero free energy: 

T = —,{2h 2 -g 2 ) + ... 
m 6 

Substituting this into ( |2.5| ) gives the effective superpotential 

Weff = NSlog(S/m s ) - 2Tvir S + 3(g 2 - 2h 2 )N^r + ... 

m 6 

Finally, extremizing it with respect to S we obtain the value of the effective superpotential 
in the vacuum: 

W eff = Nm 3 (q + 3(g 2 - 2h 2 )q 2 + . . . ) 

The same technique applies to any M = 1 theory that admits a large N limit. In 
particular, one can systematically compute instanton corrections to the effective super- 
potential in large class of M = 1 theories with any number of adjoint fields and generic 
tree-level superpotentials. 
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Appendix A. Large M Solution of the Two-Cut Matrix Model 

The results of perturbative expansions @ used in this section can be reproduced, in 
accordance with the observations of [f§, from the direct solution of the matrix model ( |2.1|) , 
as was shown in [|I|. In this appendix we review, for the sake of completeness, both the 
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matrix model derivation as well as the analytic form of the glueball superpotential in terms 
of elliptic functions. We take the cubic potential to be of the form: 



I I ■•-(<!>) ( -U - 1$ 3 



tr ±- $ ± - 2 - - $ ± -) =F — I 
2 K 2 J 3 K V \2) 



(A.1) 



The last line is the expansion around each of the two symmetric extrema of the potential. 
Note that we set here A = 1. 



In terms of the eigenvalues, Using (|2.13| ), we write the usual for the one matrix model 
saddle point equation (SPE) in the large M limit, in terms of the eigenvalues, 



If 1 

x 2 = 2\4- dup(u) . 

4 J x — u ' 



(A.2) 



where A = g s M is the overall 't Hooft coupling. The two-cut solution can be found in 
terms of the analytical function 



G(x) 



X 2 fX4 
+ 

x\ J x 3 



dup(u) 



1 



x — u 



1 

A 



X 



- — a/(x — Xi)(x — X2){x — X^)(x — X4) 



(A.3) 



having the large x asymptotics G(x — > 00) = 2 Sl ^f 2 and the corresponding couplings on 
each of the two intervals Sj = g s Mj, j = 1, 2, finite in the limit g s — > 0, M, Mj — > 00. 
The limits Xi, i = 1, 2, 3, 4 are defined by the large x asymptotics: 



J2 x * = ° 

i 
i 



(A.4) 



and by the normalization condition for the two intervals. The latter is given in terms of 
the elliptic integrals 



Si 



1 



2tt 



14 



dx \J (xi — x)(x — X2)(x — xz){x — 24) 

S2 = / dx y/ {x — Xi)(x — X2){x — Xs)(x — X4) 

2yr Jx 3 
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(A.5) 



Let us now compute the free energy J-"o(Si, S 2 ) = logZ. From the eigenvalue represen- 
tation of the matrix model we obtain the derivative of the free energy, amounting to the 
removal of the eigenvalue at the edge of each cut: 

ds 1 J 7 o(Si, S 2 ) = 97 1 (Fo(Si, S 2 ) - Fo(Si - g s , S 2 )) = -^ W M + j^t lo §( x i ~ 



and a similar expression for ds 2 J~'o(Si, S 2 ). In terms of the eigenvalue density this gives: 

\ 2 d Sl MSi,s 2 ) = 

<£ + (v \J [x\ — x)[x — x 2 )(x — xs)(x — X4) log(xi — x)dx 

.J x\ Jx 3 . 



W(r, l) + - 



= WM + Si 



■''2 



+ 



■t-'i 



■r:i 



\J (x~i — x){x — x 2 )(x — xz){x — X4) log(x4 — x)dx 



By expanding the contour of integration we pick up the contribution on the logarithmic 
cut (apart from singularities at x = 00 which we have to subtract in the matrix model 
framework). This gives: 



where Aq 



ds 1 ^o(Si : S 2 ) = W(xi) + 111 + subtractions for Ao — > 00 
ds 2 ^o(^i^ $2) = W(x4) + n 2 + subtractions for A — > 00 

00 is a cut-off and 



(A.6) 



1 



Xl 



K J -A 



yjx — xi)(x — x 2 ){x — xz){x — x^jdx 



n 9 = -■ 



TV 





Ao 



\J {x — xi)(x — x 2 ){x — Xsjjx — x^jdx 



(A.7) 
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are the dual periods. Formulas of this type appeared in |31[ , see also |jj . In they follow 
from the analysis of the Calabi-Yau geometry with flux. Using (|A.5|) , (|A.6| ) and (|A.7| ) and 
( |A.4|) one finds the small S± , ^2 expansion for the free energy itself (|3J]) from . 

Let us note that the branch points are not necessarily placed on the real axis. For 
a general complex g s , they will choose their positions according to the steepest decent 
in the eigenvalue integral. For a real g s the stable cut will be on the real axis, whereas 
the unstable cut will cross the real axis, having the complex conjugated branch points. 
The situation when all branch points are on the real axis corresponds to the analytical 
continuation in the (originally positive) variables: Si > 0, S 2 < 0. 
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A.l. Symmetric Filling of Two Intervals 

Let us consider the case of the symmetric filling of two intervals x G (5, a) and x G 
(—a, —b). It corresponds to the "unphysical" filling parameters \g s M = S = Si = —S2 > 
0, but nevertheless it will reproduce the corresponding particular case of planar graph 
expansion considered in the previous section. One can say that the two intervals are filled 
by M/2 eigenvalues and M/2 "holes", respectively. As discussed in section 3.1, this case 
describes the SU (2) Seiberg-Witten solution. We will also see yet another way of obtaining 
c = 1 noncritical string at a self-dual radius from matrix models, when the endpoints of 
the cuts approach each other. 

The function G(x), having the large x asymptotics G(x — * 00) = 0(1/ x 2 ), can be 
represented as 

A G(x) = x 2 - 1 - ^J(x 2 -a 2 )(x 2 -b 2 ) 

The function G(x), having the large x asymptotics G(x — > 00) = 0(l/x 2 ), can be 
represented as 

A G(x) = x 2 - ~ - ^(x 2 -a 2 )(x 2 -b 2 ) 

The large x asymptotics fixes one relation between a and b: a 2 + b 2 = and the normal- 
ization of the density f£ j^c-\/ (a 2 — b 2 )(x 2 — b 2 ) = 1 gives the relation (using |3!|, 217.270 
and 361.01): 

\S = — [ dx^{a 2 -x 2 ){x 2 -b 2 ) 

3 Jb (A.8) 

= — [(2 - m)E - 2(1 - m)K] 
6n 

where where K(m) and E(m) are the elliptic integrals of the I-st and Il-nd kind, a = ^/ 4 ^ 2m 
and the elliptic nome is m = 1 — b 2 /a 2 . The derivative of the free energy (|A.6|) can be 
calculated by the deformation of the contour to the dual period correspoding to the interval 
(—6, b) as the complete elliptic integral 

d s F(S, — S) = y / ^{x 2 -a 2 ){x 2 -b 2 )dx (A.9) 
A J-b 

However, the simplest quantity to calculate is actually the second derivative of the free 
energy, which is to be identified with the r-parameter of the SW curve. The latter can be 



beware of a mistake there: g 
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seen already in the form of (|A.9|) . Indeed, by writing (x 2 — a 2 )(x 2 — b 2 ) = (x 2 — |) 2 — A 4 , 
where A 4 = tkpt — m, we obtain 

16(2 — m) A ' 

d 2 MS, S) = &mds ^ ~ S) = 4K (1 - m)/K(m) = 4r (A.10) 

We found the explicite elliptic parametrization of the free energy: it is parametrized by m 
which can be expressed through S by means of ( |A.8| ). 

Expanding flAj) and ( |AloD in powers of A 4 = 16 ^" m p , we get 

S = ( A) V 6( jLf + 140( ^)- + 4620(^) 16 + . . . (A.11) 

i( , / i2/ „x 5 44 269 a8 1939 Al2 922253 Alfi \ 
r = -- (2 log(A 2 /8) + ^A 4 + ^ A 8 + ^A 12 + —^—A G + ...). (A.12) 

The last is precisely the instanton expansion of the the SW coupling constant, see e.g. 
||33j| . It is not surprizing since the numerator and denominator of ( A. 10 ) coincide (up 



to the same factor) with ood and u from the formula (2.2) of [33|, whose ratio defines r 
of course. Restoring the modulus u and rescaling A 2 — > 2A 2 , one can write this result 
in conventionsi of ref. |34|| . which also agree with our conventions used in section 3.1. 
Specifically, one finds (up to an overall numerical factor): 

, / . o i . ^ 5 A 4 269 A 8 1939 A 12 922253 A 16 . 
r = log(AV4u) + -- + —- + —— + —— + ... (A.13) 



This is in agreement with ( |3.8j ), as follows directly from the identification (|3.7| ) of the 
u- variable with A 2 /4. 

It is not surprizing that inverting the series for S plugging it into ( |A.10| ) and expanding 
in S we obtain 

1 42520 

d 2 s F (S,-S) =21og^ + 685+ 15005 2 + — - — S 3 + 0(S 4 ) (A.14) 

which coincides in the particular case S = Si = —S% with the expansion from quoted 
in section 3.1 (we put the dimensionful coupling A = 1). 

Using ( |A.8| ) and ( |A.10| ) we could also expand T itself in terms of the variable q = 
e~ 7TT = e - 7rK ( 1 - m )/ K ( m ) ? which will be the instanton expansion for the corresponding 
M = 1 SYM theory with the U(2) — > U(l) x U(l) symmetry breaking cubic tree potential, 
according to the recipe of 0. 



See also the footnote on page 3 in 
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A. 2. c = 1 critical regime 

In the context of the cubic potential matrix model, there are two distinct ways of 



getting a c = 1 non-critical string: As noted in |35|] c = 1 at self-dual radius is equivalent 
to topological B-model on the deformed conifold, which in turn has been shown to be 
equivalent to matrix model with quadratic potential [1]] . Thus in the theory we are dealing 
with, if we zoom to the region near the critical points of the potential we obtain a c = 1 
system at self-dual radius. However, there is another way of obtaining c = 1 as well: We 
can consider the limit where the two ends of the cuts touch each other, which again leads 
to a conifold geometry but now the vanishing cycle is "magnetic" relative to the original 
"electric" cycle of the matrix model. 

Let us now look at this regime in more detail. This corresponds to b — > 0, when 
mi = 1 — m —>■ and the two cuts in F(x) merge into one. From (|A.8|) we obtain in this 
limit 

A5 ~ — mi log(l/rai) + OOmf) 

D7T 57T 



6ttA 

is 

12 



which shows that this transition happens at S c = . The free energy ( [A.lOp in this limit 

w--*)- spkr + ° (mi) 

has, as a function of 5 = X(S C — 5), the typical behaviour of the c = 1 matter coupled to 
the 2D gravity 

MS) ~ Fo(Sc) ^ A 2 



log^ 

Appendix B. Massive Vacua of J\f = 1* Theory 

In this appendix we discuss matrix perturbation theory for non-trivial massive vacua 
of M = 1* theory, corresponding to higher spin representations of 577(2). As we shall see, 
there are some novelties in perturbation theory, which make these vacua conceptually sim- 
ilar to multi-cut matrix models. In fact, the vacua we are going to discuss also correspond 



to multi-cut matrix model JHJ . In both cases one finds (partial) gauge symmetry breaking 
which leads to new fermionic ghost degrees of freedom. 

In order to describe this more specifically, let us rewrite the tree-level superpotential 



( 3. 11|) in Af = 1* theory in the following form: 



W tree = tr $ 2 ]$ 3 + ^ 
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Supersymmetric vacua of the gauge theory correspond to the critical points of this super- 
potential. Thus, extremizing (|B.1|) we find 

[$i,$ 2 ] = 2z$ 3 (B.2) 

plus two similar equations obtained by permutation of indices 1, 2, 3. One obvious solution 
corresponds to $i = 0. However, there are also some non-trivial solutions, corresponding 
to p-dimensional representations of SU{2). In fact, suppose we start with a U(N) gauge 
theory, with N = pn. Then, we can take n copies of such p-dimensional representations. 
This leads to a partial breaking of gauge symmetry, 

U(pn) -> U(n) (B.3) 

Note that the rank of the gauge group has been reduced in this case due to the fact 
that the irreducible representation we have taken for vacuum configurations are not one 



dimensional. The exact effective superpotential for all values of p is known [30,37, 10 1 , and 
can be written in terms of the Eisenstein series E^Cr), 

Nv 2 

W eS = -^-E 2 (r) (B.4) 

very much like the superpotential in the for trivial vacuum, p = 1. The only novelty here 
is the relation between r and the bare coupling constant, 

r = p(pr + k)/N 

In the effective field theory, this relation is set by the tree-level term and the one-loop 
anomaly term in the superpotential. The functional dependence on r, on the other hand, 
is determined by matrix perturbative expansion JF (around the corresponding vacuum). 
Since for all values of p we have the same functional dependece on r - - given by the 
Eisenstein series — we conclude that Tq should be the same for all vacua, i.e. for all 
values of p: 

Fo = -S 3 + 7 -S 4 + ... (B.5) 

In order to reproduce this result directly by perturbative techniques in matrix model, 
we have to expand the superpotential (|B.1|) near a vacuum: 

+ $i^y + $ 2 , $i^Z + $ 3 
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where X, Y, and Z solve ( |B.2|) : 



[X, Y] = 2iZ, etc. (B.6) 
Substituting this into ( [B.1J ) we find: 

W tree = tr $ 2 ]$ 3 + J2 ^ + + ir[$ 3 , $1] + iZ[$ u $ 2 ]^ (B.7) 

Let us consider a specific case, corresponding to p = 2. In this case, we have the 
following gauge symmetry breaking pattern: 

U(2M) -> C/(M) (B.8) 

Hence, it is convenient to write all the matrix variables in terms of M x M blocks. Specif- 
ically, we take (it is easy to check that this is indeed a solution to ( p.2| )): 

*=.'? 1 ), Y = i lQ - 11 



i o j ' v.i / v o -l 

and for each hermitian matrix $j we introduce the notation 



2\D-iF A+ -A- 

where A^, D, and F are M x M matrices. Using this decomposition for all of the three 
matrix fields we get in total 3 x 4 = 12 matrices of size M x M: 

Af, £>x, F ls Af, D 2 , F 2 , Af, D 3 , F 3 (B.10) 

However, the gauge symmetry breaking (|B.8|) suggests that 3M 2 degrees of freedom can 
be gauge fixed to zero, so that effectively we should end up only with 9 matrix fields. This 
is precisely what one finds. 

Rewriting (|B.7|) in terms of M x M matrices gives the following quadratic (mass) 
terms 

WW = tr (IT, (^ + ) 2 + \ D * + \ F * + \ ( A 3 ) 2 + D ^ + F 2 A~ Z - D 1 Az + 

(B.ll) 

+ \{D2 - Fi) 2 + i(A 2 - - F 3 ) 2 + \{A- X + D 3 ) 
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Here, the fields in the first line have non-degenerate mass matrix. However, the fields 
in the second line appear only in certain linear combinations. Hence, their orthogonal 
combinations, 

D 2 +F 1 

Az+F 3 (B.12) 

represent massless directions and can be potentially dangerous in the matrix integral. In 
fact, these are simply the usual Goldstone zero- modes which can be removed by gauge 
fixing. We choose the following gauge, suggested by ( |B.12| ): 

D 2 = -F 1 

A- = -F 3 (B.13) 
A~=D 3 

This eliminates three out of twelve M xM matrices. For example, if we choose to eliminate 
D2, A2 and A^, we end up with nine bosonic matrices: 

A+, D u F x , A+ F 2 , A+, A~, D 3 , F 3 (B.14) 

Next, we should introduce fermionic ghost fieldss B, C. In order to do this, we note 
that under SU(2M) gauge transformation the matrix fields $^ transform as: 

5$ ~ [$, C] 

Again, we write C in the 2x2 block form, similar to ( |B.9| ): 

c -Uc D % F Cd -c a Cf ) < b - 15 ) 



Applying the gauge transformation to (|B.13| ) and using the standard Faddeev-Popov 



method, one finds the action for the ghost fields B ai C a , where we introduced a new 
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index notation a = A, D, F. Straighforward, but slightly technical calculation gives: 
W ghost = tr (8iB A C A - 4iB D C D - 4iB F C F + 

\2iC A {D x - F a ) + 2i(D 1 - F 2 )C A + +C D (-A+ - iD 3 ) + C F (-A+ - iF 3 ) 

+ 

2C A (W 3 - A+) + 2(A+ + iD 3 )C A + C D (iF 2 - iA 3 ) + C F (-A+ + iF x )+ 



+ \b a 



+ (A+ - iD 3 )C D + (A+ - iF 3 )C F 

+ \b d 

+ (iF 2 - iA~)C D + {iF x + A+)C F 

+ -B F 
2 



+ 



(B.17) 



(B.18) 



2C A (iF 3 - At) + 2(A+ + iF 3 )C A + C D (iF x + A+) + C F (-iD 1 - iA^)+ 

+ {iF x - A+)C D + {-iD x - iAq )C f J\ 

(B.16) 

Summarising, in the case of p = 2 we find a (9 + 6)-matrix model, that is a matrix 
model with 9 bosonic and 6 fermionic (ghost) fields, 

Bosonic : A+, D x , F x , A+ , F 2 , A+, A 3 , D 3 , F 

Fermionic : B A , Bo, B F , C A , Co, C F 

and with the following action 

Wtree = Wquadr + W^cubic + Wghost 

where the ghost action is given by flB.16|) . The quadratic terms of the bosonic action are 
given by ( |B.11| ): 

W qu adr =tr {Aff + 2Fl + 2F| + 2Dj+ 

i 

+ \d\ + \Fl + \ (A 3 -) 2 + D 1 F 2 + F 2 A 3 - - IM3 ) 
while the cubic interactions read 

W cuhic = tr (~ ([At, A+]A+ + ([F X , F 2 ] - [D x , F x ) - [D 3 , F 3 ]) A++ 

+ ([D 3 , D x ] + [F 3 , F X ])A+ + ([D 3 , F x ] + [F 2 , F 3 ]) - [At, F 3 }A 3 ~ + 
+ [D 3 , A+]A~ +i(-2F? - D X F 2 - F 2 D X )A~ + 
+ i (-2F 3 2 D 1 + F X F 3 D 3 + F X D 3 F 3 + 2F 2 D\ + F X F 3 D 3 + F X D 3 F 3 ) 



(B.19) 



Computation of the planar Feynman diagrams in this matrix model is expected to 
reproduce the perturbative expansion of the free energy (P-5|) . We will not pursue it 
further in this paper. 
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